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Abstract 

We compute the third-order corrections to the heavy quarkonium spectrum and 
production/annihilation rates due to the leading renormaHzation group running of 
the static potential. The previously known complete ©(m^af ) result for the heavy 
quarkonium ground state energy is extended to the exited states. After including 
the 0{a^g) corrections the perturbative results are in surprisingly good agreement 
with the experimental data on the masses of the excited T resonances and the 
leptonic width of the T(IS') meson. The impact of the corrections on the T sum 
rules and top quark-antiquark threshold production cross section is also discussed. 



PACS numbers: 12.38.Bx, 14.40. Gx, 14.65.Ha 



1 Introduction 



The theoretical study of nonrelativistic heavy quark-antiquark systems is among the earh- 
est apphcations of perturbative quantum chromodynamics (QCD) p. Its apphcations to 
bottomonium j2| and top-antitop jH] physics entirely rely on the first principles of QCD. 
In general perturbation theory can be applied for the analysis of these systems. Non- 
perturbative effects |4|5j are well under control for the top-antitop system and, at least 
within the sum-rule approach, also for bottomonium. This makes heavy quark-antiquark 
systems an ideal laboratory to determine fundamental parameters of QCD, such as the 
strong coupling constant and the heavy-quark masses rriq. 

The binding energy of the heavy quarkonium state and the value of its wave function 
at the origin are among the characteristics of the heavy-quarkonium system that are of 
primary phenomenological interest. The former determines the mass of the bound state 
resonance, while the latter controls its production and annihilation rates. 

Recently, the heavy quarkonium ground state energy has been computed through 
0{agmq) including the third-order correction to the Coulomb approximation [B^. The 
result has been used to extract from the T(15') meson mass. The properties of the 
excited states are more sensitive to the nonperturbative phenomena, and the correspond- 
ing perturbative estimates cannot be used, e.g., for the accurate determination of the 
heavy-quark mass by direct comparison to the meson masses. However, they have to be 
taken into account in the framework of the nonrelativistic sum rules [2] which is based on 
the concept of quark-hadron duality and keeps the nonperturbative effects under control. 
Moreover, by investigating the excited states with reliable perturbative results at hand 
one can test the effects and structure of the nonperturbative QCD vacuum. Still only a 
few states with small principal quantum numbers n and zero orbital momentum / are of 
practical interest. 

As far as the wave function at the origin is concerned a complete result is only avail- 
able through 0{a1) |7|8j . The correction has turned out to be so sizeable that 
the feasibility of an accurate perturbative analysis was challenged jU], and it appears in- 
dispensable to gain full control over the next order. Only the logarithmically enhanced 
0{al\n^as) [TniTT] and ^(a^lna,) fTM^ (for QED, see Refs. |14ll,^ll6j ) corrections are 
available so far. 

In this paper, we take the next step and calculate the nonlogarithmic third-order 
corrections to the wave function at the origin and to the spectrum of the excited heavy 
quarkonium states proportional to /?q, where /?o is the one-loop QCD beta-function. To- 
gether with the contributions already known, the new term allows to derive the complete 
result for the binding energy of the excited states. On the other hand the large-/3o terms 
usually constitute a considerable part of the corrections and can be used to estimate the 
unknown nonlogarithmic third-order contribution to the wave function. 

In the next section we present the 0{f3Qa^) corrections for the states with principle 
quantum number n = 1, 2, 3 and angular momentum 1 = 0. In Section |21 we generalize 
the complete C(mgQ!f) result for the ground state energy [0] to the excited states. In 
Section 13 we discuss the impact of the corrections on the phenomenology of the bb and ti 
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systems. Our summary is presented in Sectional 



2 Heavy quarkonium parameters to 0{pQa^ 



In the framework of nonrelativistic effective theory |17ll8ll9l2Uj the corrections to the 
heavy quarkonium parameters are obtained by evaluating the corrections to the Green 
function of the effective Schrodinger equation j^. The (3^ part of the third-order con- 
tribution results from the leading renormalization group running of the static potential 
which enters the corresponding effective Hamiltonian and is given by (see also Ref. 



Vc{r) 
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where L,. = ln{e'~"^ fir) , je = 0.577216... is Euler's constant, ({z) is Riemann's zeta- 
function with value C(3) = 1.202057..., Cp = (N^ - l)/(2A/'c) and Ca = for the 
SU{Nc) gauge group. Furthermore, we have as = if not stated otherwise. The 

coefficients (i = 1,2) and j3i (i = 0, 1,2) are given in Appendix 1X1 For the three- loop 
coefficient as only Fade estimates are available so far [23 • In the order of interest one 
has to consider single iterations of the /3q term, double iterations of the Pq and /?o term 
and triple iterations of the first-order corrections proportional to f3o- For the practical 
computation we use the method elaborated in Refs. |24|7|25j . In this way we obtain the 
corrections to the energy levels and wave function at the origin in the form of multiple 
harmonic sums. For general n the result is rather cumbersome. For a specific n, however, 
the summation can be performed analytically. Below we present our result for n = 1, 2, 3 
and / = which is sufficient for the phenomenological applications. For vanishing angular 
momentum we can write the perturbative part of the energy level with principal quantum 
number n as 



+ 6Ei'^ + 5^(2) + 5^(3) + 
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where SE^^ stands for corrections of order a^. The leading order Coulomb energy is given 
by 
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For the 0(/3Qa^) term we obtain 
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where L„ = ln{nfi/{CFas{fJ')mq)) and C(5) = 1.036927. . .. Note that the n = 1 result has 
aheady been known (26|27|6j. The perturbative expansion for the wave function can be 
written as follows 



where 



(5) 
(6) 



is the leading order Coulomb value. Our result for the (9(/?Qa^) term reads 
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3 Exited states spectrum to 0{mqal) 

The heavy quarkonium spectrum up to 0{mqa'^) has been derived in Refs. |28|7|8j . For 
convenience of the reader the expressions for SE^^^ and 5E^^^ are listed in Appendix |Bl 
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At 0(mqag) it is convenient to split SE^j^^ into two parts: one corresponding to vanishing 
beta-function and one proportional to the coefficients of beta-function: 



+ 5E^) 



The contribution 



/3(as)=0 



has been evaluated in Ref. [21] ■ For completeness we include 



the corresponding expressions in Appendix B. In Ref. P the quantity 5E'^'> has been 

computed for n = 1. Below we extend it to the excited states. Following Ref. [H] we 
divide 5^?) 



SEi'-) 



into four pieces 
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The first two terms of the above equation are related to the running of the lower-order 
potentials. The contribution 5E^^ ^ is due to the three-loop running of the static 
potential, Eq. (|T)). It reads 
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where P„ = \E'i(n) + 7£;, '^n{z) = d"ln(r(2;))/d2;" and T{z) is the Euler's gamma-function, 
in Eq. is included in Eq. (jH). 

is due to the one-loop running of the power suppressed 



The term 6^^} En 



c.r. 
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B.r. 



terms in the NNLO^ effective Hamiltonian (see, e.g., Ref. [ll]), which we denote as the 
"Breit potential" . For this contribution we obtain 
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where S is the spin quantum number. 



^LO, NLO, . . . stand for the leading order, next-to-leading order, etc. 
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The remaining two contributions of Eq. Q are related to the iteration of lower-order 



potentials. The contribution 6E^^ 
loop running of the static potentia 

3 /■ 
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where ^Ls 

'^0 



c.i 



contributes to Eq. (jl)). 



The last contribution 5Ei^^ 



the one-loop running of the static potential. It reads 
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After summing up the four contributions according to Eq. © we obtain our final result 
for the 0{al) corrections to the energy levels involving coefficients of the beta function. 
For n = 1, 2, and 3 they read 
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The equation with n = 1 agrees with the result of Ref. 0. The Eqs. (HD), (|^. (pUj) 
and provide the complete result for the energy levels up to O(mqa^). We should 
note that, although we only present analytical results for the first three principle quantum 
numbers, there is no principle problem to obtain expressions for higher excited states, too. 
However, from the phenomenological point of view they are far less important, and thus 
we refrain from listing them explicitly. 

It is instructive to evaluate the energy levels in numerical form: 

6E?^ 



6E. 



(3) 



a. 



70.590|„,=4 
56.732|„,=5 

84.634|„,=4 
62.164|„,=5 



15.297 In a, + 0.00103+ I V^'ITa i 

\ 26.654|„,=5 



+ 8.6471na, + 0.00103 



f 34.2291, 
^ 26.654], 

67.337|„,=4 
52.434U,=5 



03 



101.69U,= 
72.368L,= 



+ 6.305 Ina, + 0.001 aa + 



ni=4 \ 




ni=5 J 


Pi. 



(15) 



where as = as{ps/n) and /i = ps/n with ps = CFOis{ps)'mq and we put S = 1 which 
corresponds to the spin-triplet state. The recent analysis of the spin-dependent contri- 
bution to the spectrum, which is responsible for the hyperfine splitting, can be found in 
Refs. |29|30j . In Eq. (fT3j) we have separated the contributions arising from 03 and /Jq. 
Using the Fade estimates j2S] we obtain 0.001 a3|„j=4 ~ 6 and 0.001 a^\ni=5 ~ 4. Thus, the 
result for the energy levels depends only marginally on the precise value of 03 provided 
the Fade estimates give the correct order of magnitude. Furthermore, one can see that 
the /?Q term contributes between 25% (n = 1) and 50% {n = 3) of the nonlogarithmic 
term. 



4 Heavy quarkonium phenomenology 

In this section we discuss some phenomenological applications of the results derived in 
the previous parts of the paper. As input values for the numerical analyses we adopt 
as{Mz) = 0.118, and mt = 5.3 GeV and rrit = 175 GeV for the quark pole masses. 
Furthermore, we use the soft scale /x^ ~ 2.10 GeV for the bottom and ps ~ 32.6 GeV for 
the top quark case. 

Excited states of bottomonium. The mass of the T{nS) meson can be decomposed 
into perturbative and nonperturbative contributions 

Mt(„5) = 2m6 + EP-*- + 5-P-E„. (16) 

The perturbative contribution up to 0{mgal) is given in the previous sections. The 
phenomenological application of the result to the T(IS') meson mass has been discussed 
in Ref. [U]. For the exited states let us consider the ratio 

En -El 

It depends on the quark mass only through the normalization scale of and does not 
suffer from renormalon contributions. Including successively higher orders one gets for 
p = Ps 

lO^xpP-*- = 1.49 (1 + 0.79nlo + 1-18 NNLO + 1-21n3L0 + • • •) , 

10^ X pf- = 1.77 (1 + 0.92NLO + 1.37nnlo + 1.55n3lo + •••), (18) 

where a^f'{iis) is extracted from its value at Mz using four-loop beta-function accom- 
panied with three-loop matching^. Though the convergence of the series is not good. 



^We use the package RunDec '31' to perform the running and matching of 
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T{2S) 


T{3S) 


10^ X p^^P 


5.95 


9.46 



Table 1: Perturbative versus experimental results for the parameter p„ as defined in 
Eq. (fTTj) . The theoretical uncertainty corresponds to as{Mz) = 0.118 ± 0.003. The 
experimental values are extracted from Ref. [HI]. For 03 we used the Fade estimate [2S| 
«3|„,=4 = 6272. 



the N^LO perturbative result is in impressive agreement with the experimental values 
Pn'^ — (^T{n5) " My (IS)) / My (IS) for n = 2 and 3 as can be seen in Tab. [TJ We would 
like to emphasize the role of the perturbative corrections necessary to bring theory and 
experiment into agreement which we will use in the following to estimate the order of 
magnitude of the nonperturbative effects. In fact the absence of a sufficiently accurate 
estimate of the nonperturbative part d'^'^'En is one of the main problems in the theory of 
heavy quarkonium. In the limit a'^uig S> ^qcd it can be investigated by the method of 
vacuum condensate expansion |4|5j . However, for bottomonium it can only be used for 
n = 1. For higher states the leading term due to the gluonic condensate grows as n^. 
It becomes unacceptably large already for n = 2 where the whole series blows up [22] • 
Even for n = 1 such an estimate suffers from large uncertainties due to the poorly known 
value of the gluonic condensate and due to a strong scale dependence. A rough numerical 
estimate is ^ 60 MeV 0. Since our perturbative result agrees very well with the 

experimental result we can conclude that 5" P i?2 should be of the same size as 5" P i?i. 
In general for bottomonium the nonperturbative corrections appear to be rather moder- 
ate and the theoretical estimates are dominated by perturbative contributions. Similar 
conclusion has been made in Ref. jHS] in a somewhat different framework. 

T(IS') leptonic width. In the nonrelativistic effective theory the leading order approxi- 
mation for the leptonic width rLO(T(l^) /+/") = Ti reads = ATTN^Qla^\t{j^{0)\^/ 
{3ml), with A^'c = 3 and Qb = —1/3. Combining the known perturbative results up to 
O^a^lnas) (see Ref. jI2I) with the 0{PQal) contribution obtained in SectionElwe obtain 
the following series 

Fi ^ (l - 1.70 asirub) - 7.98 al{mb) + . . .) 
X (l - 0.30 as - 5.19 In a, + 17.2 

-14.4a^ln2a, + 0.17a^lna, -34.9a^|^3 + ...) , (19) 

where as = asips)- The contribution coming from the hard virtual momenta region |35|36j 
is separated and the corresponding strong coupling is normalized at fi = nib. Evaluating 
Eq. p9|l and retaining only the logarithmic and (3^ terms at N^LO we find 

Fi ^ rJ^O(l- 0.445NLO + 1.75NNLO-1.67N3LO' + •••), (20) 



9 



0.6 
0.4 
0.2 




4 4.5 5 
\x (GeV) 



0.6 
0.4 
0.2 




10 20 30 40 50 60 70 



90 100 
|I (GeV) 



(a) 



(b) 

Figure 1: (a) Fi normalized to Fi = ^I'^la^^asitis) ^ function of fi at LO (dotted), NLO 
(dashed), NNLO (dotted-dashed) and N^LO' (full line). The horizontal line corresponds 
to the experimental value F"^p(T(1^) ^ e+e~) = 1.31 keV For the N^LO' result, 
the band reflects the errors due to as{Mz) = 0.118 ± 0.003. (b) The analog plot for Ri 
with R, ^ R\X.-.M,sr 



where the prime indicates that the N^LO corrections are not complete. Though the 
perturbative corrections are huge, the rapid growth of the perturbative coefficients stops 
at NNLO if we assume that the (3q term sets the scale of the nonlogarithmic third- 
order contribution. In Fig. HJ^a), the width is plotted as a function of /i including the 
LO, NLO, NNLO and N^LO' approximations along with the experimental value. For 
the numerical evaluation we extract a^J^\mi,) from its value at Mz using four-loop beta- 
function accompanied with three-loop matching, a^'^^rrih) is used as starting point in 
order to evaluate a^'^^(/i) at N'^LO with the help of the {k + l)-loop beta-function. As 
one can see in Fig. ^a), the available 0{al) terms stabilize the series and significantly 
reduce the scale dependence. At the scale fi' ~ 2.7 GeV, which is close to the physically 
motivated scale fig, the N^LO' corrections vanish and at the scale n" ~ 3.1 GeV the result 
becomes independent of fi; i.e., the N^LO' curve shows a local maximum. In the whole 
range of fi between 2 GeV and 5 GeV the result for the width agrees with the experimental 
value within the error bar due to the uncertainty of the strong coupling constant. This 
may signal that the missing perturbative corrections are rather moderate. Furthermore, 
this result constitutes a significant improvement as compared to the NLL approximation 
discussed in Ref. [Tf] . 

For a definite conclusion, however, one has to wait until the third-order corrections 
are completed. The potentially most important part to be computed is the ultrasoft 
contribution which includes Q;s(/i) normalized at relatively low ultrasoft scale fius ~ alrriq. 
Currently only a partial result for this contribution exists |38j. 

T sum rules. The nonrelativistic T sum rules j2] operate with the high moments of 
the spectral density with n ~ l/«s, which are saturated by the nonrelativistic near- 
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threshold region. The experimental input is given by the masses and leptonic width 
of the T resonances which are known with high accuracy. On the theoretical side the 
nonperturbative effects are well under control. This makes the T sum rules one of the 
most accurate sources for the bottom quark mass value. The complete perturbative 
analysis has been performed up to NNLO |24p7pl6 39,27] . The extension to N^LO is a 
challenging problem. 

The theoretical value of the high moments is saturated by the contribution of a few 
lowest heavy quarkonium states and the corrections to the moments are dominated by 
the corrections to their masses and wave functions at the origin. To estimate the size 
of the N^LO corrections we include the 0{mqal) result for the energy levels and the 
partial 0{al) result for the wave function at the origin which includes all the logarithmic 
term J2] and the (5^ terms obtained in Section El We perform the analysis along the 
lines described in Ref. j7] using n = fig. For n > 20 the corrections to the moments are 
dominated by the one to the ground state energy and we recover the result of Ref. jB] 
for the bottom quark mass. For lower moments, which provide better balance between 
theoretical end experimental uncertainties [7], the situation changes drastically as the 
corrections to the wave function at the origin begin to play an important role. For n = 4 
the negative third-order contribution to the wave function completely cancels the effect 
of the third-order correction to the binding energy, and the correction to the pole mass 
rub almost vanishes. The pole mass can be converted into the MS mass rhb{rhb) which is 
widely believed to have much better perturbative properties. If we correlate the series so 
that the A;*'^-order correction to the sum rules goes along with the fc-loop mass relation, 
which is natural for low moments, we obtain as an effect of the third-order corrections 
6Thb{mb)^3i^o PS —100 MeV. We take this variation as an estimate for the size of the 
N^LO corrections within the T sum-rule approach. It is interesting to note that the 
N^LO correction to fhb{ffib) is negative at the soft normalization scale in contrast to the 
series obtained from the ground state energy analysis 0- 

Top quark-antiquark threshold production. The nonperturbative effects in the case 
of the top quark are negligible. However, due to the relatively large top quark width, 
Tti its effect has to be taken into account properly [Hj since the Coulomb-like resonances 
below threshold are smeared out. Actually, the cross section only shows a small bump 
which is essentially the remnant of the ground state pole. The higher poles and continuum, 
however, affect the position of the resonance peak and move it to higher energy. The value 
of the normalized cross section R = a{e~^e~ — > tt)/a{e^e~ Z^'*'/^") at the resonance 
energy is dominated by the contribution from the would-be toponium ground state which 
in the leading approximation reads R\'^ = GvrA^cQt iV'f (0)P/ [nT't^t), where Qt = 2/3. 
The analog to Eq. (fT^ reads 



Ri ^ R^^ (l - 1.70 asimt) - 7.89 a^(mt) + . . .) 




13.8 In^ a, + 2.06 In a, - 27.2 a^|^3 + . . .) , 



(21) 
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with as = Numerically we find 

i?i ^ i?5^o^^_o.243NLO + 0.435NNLO- 0.268N3LO' + •••)• (22) 

The new third-order corrections proportional to /5g amount to approximately —7% of the 
LO approximation at the soft scale which is the same order of magnitude as the 0{al) 
linear logarithmic term. The available N^LO terms improve the stability of the result 
with respect to the scale variation as can be seen in Fig. ^^b). The absence of a rapid 
growth of the coefficients along with the alternating-sign character of the series and the 
weak scale dependence suggest that the missing perturbative corrections are moderate and 
most likely are in the few-percent range. It is interesting to note that the perturbative 
contributions of different orders, which are relatively large when taken separately, cancel 
in the sum to give only a few percent variation of the leading order result. 

5 Summary 

In this paper the important class of the third-order corrections to the heavy quarkonium 
parameters proportional to /5q has been obtained. The complete result for the exited 
states spectrum to 0{mga^) is derived. The perturbative results are in surprisingly good 
agreement with the T(2S') and T{3S) meson masses and the leptonic width of the T(IS') 
meson. Thus the nonperturbative effects in bottomonium seem to be rather moderate 
and the theoretical results are dominated by the perturbative contributions. A failure of 
early low-order perturbative analysis to describe the T system is due to large perturbative 
corrections to the Coulomb approximation. On the basis of our results the magnitude of 
the N^LO corrections to the T sum rules and top quark-antiquark threshold production 
cross section is estimated. The available N^LO corrections which include all logarith- 
mic terms and the nonlogarithmic (3^ contribution stabilize the perturbative series for 
the production/annihilation rates that makes us more optimistic about possible accurate 
perturbative description of these quantities. 
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A Static potential and beta-function 



For convenience of the reader we list in this appendix the result for the coefficients of the 
static potential (see |41f42p43j and references therein) 



31 



^2 



20 



Tpui 



4343 



22 



+ 47r^ - — + — C(3) 
162 4 3 ^ 



cl 



1798 56 
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55 
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16C(3) 



CpTpni + ( Y^F^i 



(23) 



and the beta-function 
1 



^2 



4 



34. 



16 
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44 



20 

y 



1415 



205 



158 



9 



27 



where Tp = 1/2 and ri; is the number of the light quark flavours. 



(24) 



B Results for 5E^^ 

In this appendix we collect the known results for the perturbative corrections to the heavy 
quarkonium spectrum. The first and the second order corrections read ;28^7j5] 
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(25) 



E? ( — 

TT 



l2PlLl + (3ai/?o + 4A + (-8 + 24P„+i) ^l) L„ + ^ + | 



+ (-1 + Wn+i) ai/3o +{ — + ^- 



- 16^2(^) - An^3{n)jP^o + 4P„+i/3i 
2 ^2 11 2 



H CaCf 

n 



The result for 5E^^ 



/3(as)=0 



,n IQn? 3n 
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(26) 
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.oo o 
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ClC, 



n 



139 
47 2 



7 41 
41n2 + -(7i5-lnn + ^i(n + l)) + — 

D 
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+ 
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107 



108 12n 6 
79 7 
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7 



, -ln2 + -(7s-lnn + ^i(n + l)) + 3L„ 
18 on 3 3 
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S{S+1) 



CI 

n 



32 

+ 21n2 + (l-ln2)5(5 + l) 

15 



n 



^mCACFTFTM^ 

36n 



10 



9 18n 27 



5(5+1) 



CITfUi , 2 



n 



(27) 



where 
values 



\n{CF<ys) and stands for the QCD Bethe logarithms with the numerical 



Lf = -81.5379^ 



-37.6710 



-22.4818. 



(2^ 



The terms proportional to L^^ have been computed for the first time in Ref. 
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